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We consider the slow flow of a viscous incompressible liquid in a channel of constant but arbitrary 
cross section shape, driven by non-uniform suction or injection through the porous channel walls. 
A similarity transformation reduces the Navier-Stokes equations to a set of coupled equations for 
the velocity potential in two dimensions. When the channel aspect ratio and Reynolds number 
are both small, the problem reduces to solving the biharmonic equation with constant forcing in 
two dimensions. With the relevant boundary conditions, determining the velocity field in a porous 
channels is thus equivalent to solving for the vertical displacement of a simply suspended thin plate 
under uniform load. This allows us to provide analytic solutions for flow in porous channels whose 
cross-section is e.g. a rectangle or an equilateral triangle, and provides a general framework for the 
extension of Berman flow (Journal of Applied Physics 24(9), p. 1232, 1953) to three dimensions. 



INTRODUCTION 



Channel flows - liquid flows confined within a closed conduit with no free surfaces - are ubiquitous. In animals 
and plants [f| they serve as the building blocks of vascular systems, distributing energy to where it is needed and 
allowing distal parts of the organism to communicate. When constructed by humans, one of the major functions of 
channels is to transport liquids or gasses, e.g. water (irrigation and urban water systems) and energy (oil or natural 
gas) from sites of production to the consumer or industry. 

In some cases, the channels have solid walls which are impermeable to the liquid flowing inside. In other cases, the 
channels have porous walls which allow the liquid to flow across the wall and thus modify the axial flow. Both are 
important. The first class of flow has been studied in great detail, and analytical solutions are known in a few, but 
important, cases The latter class has received much less attention, although it is equally important. The effect 
of porous walls is especially important in the study of biological flows due to the presence of permeable cell walls @ 
and in industrial filtration applications [Io| . 

Some analytic solutions of the flow in porous walled channels are known, primarily due to a similarity technique 
first used in this context by Berman [3j. Berman's method is closely related to those commonly used in boundary 
layer theory [l3j and allows for the solution of steady flows in geometries with symmetries which makes the problem 
two-dimensional. By demanding that the solution be of similarity form the Navier-Stokes equation is reduced to a 
single non-linear third-order differential equation for the velocity potential in one space dimension. The flow between 
parallel plates Q and in a cylindrical [l7j and annular tube 0] have been analyzed in this way. Time dependent flows, 
flows at high Reynolds numbers and questions of uniqueness and stability of these flows have since been address by 
a large number of workers, see e.g. [Ej]. 

In this paper, we extend Berman's method to three dimensional similarity flows, and derive a set of equations for 
the velocity potential which are valid in channels of arbitrary cross section shape. At low Reynolds numbers, and 
when the channel is very long compared to its characteristic transverse dimension, the Navier-Stokes equation reduces 
to a single partial differential equation for a velocity potential in two space dimensions; the inhomogenous biharmonic 
equation with constant forcing. This equation, which is derived in Sec.„ has been widely studied in the literature as is 
describes the transverse displacement of a simply suspended thin elastic plate under uniform load. In Sec.„ we provide 
analytic solutions to four cases of porous channel flows in geometries where the solution of the corresponding elastic 
problem is known: Flow in a cylindrical tube, between parallel plates, in a triangular channel, and in a rectangular 
channel. 



FLOW IN CHANNELS WITH POROUS WALLS 



We consider a long, straight channel parallel to the z-direction, and assume that it is translational invariant along 
this axis with an arbitrary, but constant, cross section ft as shown in Fig.QJa-b). The channel has length L, perimeter 
S, cross section area A, and volume V = LA. The coordinates in the transverse xy-plane are denoted r±_ = (x, y), so 
that the full coordinates are written as r = r± + e z z and likewise for the gradient operator V, Laplace operator V 2 , 
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FIG. 1. Flow in channels with porous walls, (a) Schematic view of the channel cross section geometry, (b) Schematic view of 
the axial flow geometry. See Sec. „ for further details, (c-f) Sketch of geometries considered in Sec.„ In (a) and (c-f) the z-axis 
is pointing into the page. 



and velocity field u 

r = r± + e z z where r±_ = (x, y) 
V = V± + e z 8 z V± = (d x ,d v ) 

V 2 = Vi + d 2 Vi = d 2 x + d 2 

u = u±+e z u z u ± = (u x ,u y ) 

where we have used the short hand notation d x f = df/(dx). 

We consider the case of incompressible Newtonian fluids of viscosity r\ and density p (in the laminar regime) which 
are governed by the Navier-Stokes equation 

p (d t u + (u ■ V)«) = -Vp + t]\7 2 u (1) 
V-m = (2) 

where p is the pressure. Writing the velocity field u in decomposed form as u = u±_ + e z u z Eqns. ([I} and @ are 

p (d t u z + (u± ■ Vi + u z d z )u z ) = -d zP + //(Vi + d 2 z )u z (3) 
p (d t u± + (u ± ■ V± + u z d z )u±) = -V X P + »?(Vi + d 2 z )u ± (4) 
V ± -u ± + d z u z = (5) 

We assume that the flow is driven by a prescribed injection or suction of fluid through the porous wall which leads 
to a normal flow velocity component at the channel wall u w (r) = w(r±)f(z) of typical magnitude u$. The boundary 
conditions thus require that the tangential velocity component vanishes on the channel wall and that the normal 
velocity component is u w [r) 

u z (r) = 0, re On, (6) 

t-u±(r) = 0, re dn, (7) 

n-u±(r) = u w (r) = w(r±)f(z), r e dn, (8) 



where n is a outward pointing normal unit vector and t is a tangent unit vector to the boundary <9f2 in the r± plane. 
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Non-dimensionalization 

To cast Eqns. (U)-© into a simpler form we non-dimensionalize the equations using the characteristic wall flow 
velocity uq, channel length L, and transverse dimension a = A/S: 

2, 



z — Lz r± — ar t = — t u± = uqu , u z = —uqu? p = 5 — p 

uq a a A 



v ± = -V x vl = ±V ± 2 d z = \d z , d 2 = ±d 2 z , d t = ^d t , 

With this change of variables, Eqns. ©-© are 

^ {d v u' z + (u' ± ■ V' ± + u>d z ,)u> 2 ) = -d z , P > + (vl> + (f ) 2 < (9) 
^ (d t u' ± + (u' ± ■ V ± + u' z d z ,)u' ± ) = W + (|) 2 (Vi' + u' ± 

(10) 

V ± -u' ± + d z ,u' z = Q (11) 

Introducing the Reynolds number Re = based on the wall velocity uq and the aspect ratio a = , and dropping 
the primes for case of reading we finally have that 

Re {8 t u z + {u ± • V± + u z 8 z )u z ) = -d zP + (Vi + a 2 d 2 z ) u z (12) 
a 2 Re (d t u ± + (u ± ■ V_l + u z d z )u ± ) = -VxP + a 2 (V 2 ± + a 2 d 2 ) u ± (13) 
Vjl • u x + d z u z = (14) 

with the boundary conditions given in Eqns. (JH])-©- 



Similarity solutions 

The form of the boundary conditions is such that the in-plane velocity u±_ could be irrotational and proportional 
to f{z) everywhere (see Eqns. (jTj)- (JHJ) ) , while the axial velocity u z should be proportional to the total volume of liquid 
which has entered the channel at z' < z, i.e. u± cx F(z) = L z f(z') dz' + Fq. It is worthwhile to enquire if the 
differential equation permits solutions of this form, and we therefore write the velocity it as a similarity solution 

u± = -f(z)g{r x ) (15) 
u z = F(z)h(r x ) (16) 

where g(r±) and h(r±) are unknown functions of the radial position rj_ only. Using Eq. (|15[) . we find from the 
continuity equation () 14[) 

V ± .g = h. (17) 

Since u±, and therefore g, is irrotational by assumption, we may write g = where ^ is a velocity potential. 

Eq. (fl7|) then implies that h = V 2 </>. The velocity field we seek is thus of the form 

u x = -f(z)V ± <t>{r ± ) (18) 
u z = F(z)V 2 ± <P(r ± ) (19) 

Substituting Eqns. (THJ) and jTSJ into Eqns. {12J and (HSJ) we find 

Re [d t {FVx<t>) + JF ((V 2 ^) 2 - (V x ^ • V±)V ± (j>)] = -d zP + FVl + a 2 d z fVl<t> 

(20) 

a 2 Re [-d t (/V_i» + / 2 (V ± </> • Vi)Vi^ - Fd z f^j\cpV x< f>] = -V ±P - a 2 (fV 2 x V±<t> - a 2 d 2 JV x cj>) 

(21) 
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The boundary conditions for (r_i_, z) E <9£1 in Eqns. ©-([7]) are 

F(z)Vi0(r ± ) = (22) 
-(t-V±cp(r ± ,z))f(z)=0, (23) 
-(n-V ± cf>(r ± ))f(z)=w(r ± )f(z) (24) 

These boundary condition may be considerably simplified by noting that Eq. ((23)) implies that 4>(r±, z) = k is constant 
on the boundary dCl. (If d£l consists of several physically separate boundaries, k may take on different values on each 
of these). With this, Eqns. ([22|)-([24]) become 

€<P = (25) 
4> = k (26) 
-5„0 = w (27) 

where we have assumed that f(z) ^ 0, F(z) ^ 0, and use the notations d n and d 2 for first and second order normal 
derivatives. 



The case Re < 1 and a < 1 



When both the Reynolds number Re <C 1 and aspect ratio a«l are small, we find from Eqns. (|2U)) - (l2Tt that 

d zP = FVU, (28) 

V ±P = 0. (29) 

Eq. (|2"9"|) implies that the pressure is a function of z only, and we may write p — p(z). By introducing P(z) — d z p z /F, 
we find the governing equation for the velocity potential <p 

Vl0 = P (30) 

with the boundary conditions given in Eqns. (|25 |) - (|27| . 



Analogy with the theory of simply suspended plates 

The equation of motion for a thin suspended plate under a uniform transverse load q is the inhomogeneous bihar- 
monic equation 

ViW = ± (31) 

where W is the displacement at any point from the position of equilibrium and D depends only on the tension and 
mass of the plate. This equation is of the same form as Eq. (|30[) . and the amplitude of the displacement may be taken 
to represent <fi if 

P=± (32) 

It appears therefore that if a solution of the problem of a suspended plate has been obtained, a problem of viscous 
motion in a porous tube has also been solved. This analogy is exact, so long as the boundary conditions in Eq. (|25|) 
(|27p are also fulfilled in the plate problem. This occurs when the plate is simply suspended, i.e. when on the boundary 
d 2 W/dn 2 = and W = k, in which case Eqns. ([23)1 and (|2"o) are fulfilled. The final boundary condition (Eq. ([27])) 
which determines the angle of deflection dW/dn or the normal flow velocity across the membrane u w is set by the 
choice of P = q/D if f2 is simply connected or by P and k on each of the boundaries if £1 is multiply connected. 

This observation that Eqns. (|3"Tj) and (|30|) are of the same form adds to a list of viscous flow problems which may 
be solved by studying displacements and vibrations of thin plates (see e.g. Lauga et al. Q, Meleshko Rayleigh 
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Richards [12j, Taylor [14j; 
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(a) (b) 




FIG. 2. Flow in a cylindrical tube of radius h. (a) Axial flow velocity u z given in Eq. (|38|l plotted as a function of radial 
position r. (b) Radial flow velocity u r given in Eq. (|37p plotted as a function of radial position r. Note that the radial velocity 
has a minimum at r — y/2/3h ~ 0.817 ft where the value is M r (y / 2/3ft) = — a/32/27uj/ ~ —1.09 wf(z). 



FLOW SOLUTIONS 



Solutions of the inhomogeneous biharmonic equation (|30j) must satisfy the governing differential equation and 
boundary conditions characterizing each geometry. The fulfillment of the boundary conditions often presents consid- 
erable mathematical difficulties and thus in general, analytic solutions are rare. Taking advantage of the solutions 
known from plate bending theory (see e.g. Timoshenko |15| . Ventsel & Krauthammer |l6[), we are able to provide 
analytic solutions to porous channel flows in geometries where the solution of the corresponding elastic problem is 
known, such as in channels of rectangular and triangular cross section shape. In these cases the solution relies on either 
a parametrization of the boundary or a series solution. Before we consider these three dimensional cases, however, 
we illustrate the solution technique on a few two dimensional flow problems which have been solved by other means 
in the literature. 



Flow in a cylindrical tube 

Consider the flow in a cylindrical tube of radius h as sketched in Fig. Hfc). Let the normal velocity at the wall be 
u r (h) = —wf(z). Assuming rotational symmetry of the velocity field, the governing equation (|30|) is 



2„, 1 „ 1 



which has the solution 



d 4 + -d 3 r + c£ + ^d r )cj> = P (33) 

<^ = |f 4 + \( 2B i - B ^ 2 + + log r + B A . (34) 

With k — at the r = h channel wall, the boundary conditions in Eqns. (|25l) - (f!?7|) are 

-BJrdMh)) = 0, 6(h) = 0, -d r 6(h) = -w (35) 
r 

Assuming further that the radial flow component vanishes at r = (u r (Q) = 0) such that — d r 6(0) = we find for the 
for the velocity potential 6 in Eq. (|3"3"|) 

^) = (-T + T-4^h (36) 

The radial and axial velocity components can now be found from Eqns. (jT5J) - (jT^J) 

r (r 2 — 2h 2 ) w , v 

u r = -fd r 6 = -± h3 > /, (37) 

1 A(h 2 — r 2 )m 

u z - F-d r {rd r 6) = ( ' F (38) 
r 
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(a) (b) 




FIG. 3. Flow between parallel porous plates located at x — and x = h with boundary conditions u x (0) = wf and u x (h) = —wf. 
(a) Axial flow velocity u z given in Eq. (|46[) plotted as a function of vertical position x. (b) Transverse flow velocity u x Eq. (|45[) 
plotted as a function of vertical position x. 



in agreement with e.g. Aldis The two velocity components are plotted in Fig. [5] We note that the radial velocity 
has a minimum at r = y/2/3h ~ 0.817/i where the value is u r (-\/2/3h) = — y // 32/27wf ~ —1.09 wf(z). 



Flow between parallel plates 



Consider the flow between two infinite parallel plates located at x = and x — h as sketched in Fig. [IJd). Let the 
normal velocity at the walls be u x (0) — W-f(z) and u x {h) — w+f(z). The governing equation pop is 

%<f> = P (39) 

which has the solution 

p 

(j) = A + A lX + A 2 x 2 + A 3 x 3 + — x 4 . (40) 

Applying the boundary conditions in Eqns. (|25 p -(|27 p 

9^(0) =0, dl<t>{h)=Q, 0(0) = 0, -d x cj)(0) = u>_, and - d x <fth) = w+, (41) 

yields 

(w_ — w + )x 3 (w+ — W-)x A 



w-x. (42) 



h? 2h 3 

The transverse and axial velocity components can be found from Eqns. ([T8j) - ([T9"1) 

u x = = ( ^^'^ + 7 +)x3 + «,_) / (43) 

«. = F(*)flg* = 6(W -'7 3 )( ^" X) ^ (44) 
The solution for the special case w_ = — w + = w is shown in Fig. [3] In that case the velocity components are 

h 3 - 6hx 2 + 4x 3 



12{h - x)x 



h 3 

These were first by obtained Berman [3J who considered the case / = 1. 



-wf, (45) 
wF. (46) 
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(a) u z (b) u± 




FIG. 4. Flow in a porous tube whose cross-section is an equilateral triangle of height h and side length 2h/y/3. (a) Contour lines 
for the axial flow velocity field u z (x, y). The contour lines are shown in steps of 10% of the maximum value, (b) Stream(thick)- 
and contour(thin) lines for the transverse velocity field u± — (u x ,u y ). The contour lines show the velocity magnitude \u±\ in 
steps of 10% of the maximum value, (c) A plot of u z (x,0) along the centerline parallel to the x-axis. (d) A plot of u x (x, 0) 
along the centerline parallel to the a;-axis. 



Flow in an equilateral triangle 



Consider the flow of liquid in a porous channel whose cross-section is an equilateral triangle of height h and side 
length 2/i/v3, as shown in Fig. HJe). Let the normal wall velocity be given by u w — w(x,y)f(z), and let the mean 
value of w(x,y) when averaged over the channel wall be wq. The governing equation (|30j) is 

{dl + d^cp = P, (47) 

while the boundary conditions in Eqns. (P5|) - (P7|) are 

dfy = 0, = 0, and - d n (j> = w. (48) 

The solution of Eq. (|47|) with the boundary conditions d^<j) — and <fi = is 

= e4 ( x3 " 3y2x ~ h{x2 + y2) + Ti h *) {t h2 ~ x2 ~ y2 ) (49) 

The magnitude of the wall normal flow velocity u w = w(x, y)f(z) is the same on each of the three boundaries. At the 
x = —h/3 wall we find 

P(h 2 -3y 2 ) 2 

u x (-h/3, y) = -fWlx, y)\ x= _ h/3 = ~jQ2h f - (50) 
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(a) u z £ = h (b) u ± £ = h 

I £/2 x 




FIG. 5. Flow in a rectangular porous tube with sides of length h and I. (a) Contour lines for the axial flow velocity field u z (x, y) 
for the case h — I. The contour lines are shown in steps of 10% of the maximum value, (b) Stream(thick)- and contour(thin) 
lines for the transverse velocity field ui for h = I. The contour lines show the velocity magnitude \u±\ in steps of 10% of the 
maximum value, (c) A plot of u z (x, 0) (solid line) along the centerline parallel to the x-axis for £/h — 1000. Dashed line shows 
u z (x, 0) for parallel plates, c.f. Eq. (|46[) . (d) A plot of u x (x, 0) (solid line) along the centerline parallel to the :r-axis £/h = 1000. 
Dashed line shows u x (x,0) for parallel plates, c.f. Eq. (|45[) . In (a) and (b), the first 100 terms in Eqns. ()55[) - (|57|) were used. 



The mean velocity (when averaged over the x = —h/3 boundary) is —h 3 Pf/360 so by choosing P = —360wo/(h 3 ), 
the average inflow velocity becomes wof. Using P = —360wo/(h 3 ) in Eq. (|4^)l . the velocity components can be found 
from Eqns. (fT5]) - (fTT))) recalling that u x = —fd x (j), u y — —fd y <f) and u z — F(d% + 9y)(p 

5(l8(6hx + 9x 2 -2h 2 )y 2 + 8ly' l -x{8h+l5x){2h-3x) 2 ) , x 

u x = ^ -o/ (51) 

5(/i + 3a;)y(8/i 2 -6ax-9(a; 2 +3 2 / 2 )) , x 
—n-,.[ (52) 



6h 4 

10{h + 3i) ({2h - 3x) 2 - 27y 2 ) 



3/i 4 

The velocity field is shown in Fig. [4] 



w F (53) 



Flow in an rectangular channel 



Consider the flow in a porous channel whose cross-section is a rectangle with sides of length h and £ centered at 
(x,y) = (0,h/2) as shown in Fig. [TJd). Let the normal wall velocity be given by u w = w(x,y)f(z), and let the mean 
value of w(x,y) when averaged over the channel walls be wq. 

The equations of motion and boundary conditions are once again given by Eqns. (|47[) - (|4"5)l . As we need to evaluate 
the first and second order x and y-derivatives to determine the velocity field, Levy's solution of the corresponding 
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plate bending problem is convenient 

, 4P/i 4 



— I 1 - A n cosh — h B n y sinh — — I sin ■ 



7r 5 ' n 5 

71=1,3,... 



(54) 



where a„ = nn£/(2h), A„ = (a„tanha„ + 2)/(2 cosha„), and B n = a n /(£ cosha„). From Eqns. (fMf , (fT£|). and (TlT))) 
we may calculate the velocity components 



4Pft 3 



y J_ 

7r 5 n 5 



n=l,3. 



/ 51 ~J ( 1 _ ^" cosl1 + sinh 1 cos ■ TWI ' 



n=l,3. 



2a n . 2a„y / . 2a n y 2a„ 2a„y 
- — g- Ai Sinn — ^ h O n I smh — h —g-y cosh — — 



* 3 „=fe,... " ;i v 



cosh a r . 



sin ■ 



To determine P such what the average inflow velocity is wof, we solve for P in 

1/2 



- Uy(-£/2,x, z) dx + - / Uj.(0,y,z) dy = w f 
n Jo 1 J —1/2 

In the following we determine P for the special cases £ = h and £^> h. 



(55) 

(56) 
(57) 

(58) 



The case t = h 

To determine the velocity field in the case £ = h we note that the sums in Eqns. (|55[) - (|57p converge very rapidly, 
and that the flow profile u x (0, y) at the x = wall is well approximately by the first term in Eq. (|55[) 



u x (0,y,z) 
From Eq. ([55} we therefore find 



APh 3 ( 2a n y 2a„y\ r 

I 1 - An cosh — — + B n yswh ) f 



u x (0,y,z) dy • 



8h 3 P 



-1/2 



£tt\- 



i-K 



I** \ SGCh 2h (M2 + co S h- )-3h^h T ]f = wof, 



which determine P as a function of the average normal flow velocity wq 



P = - 



1 + COSh 7T 



w 



16 7r(2 + cosh7r) — 3 sinhw h 3 



-29.90^ 
h 6 



(59) 



(60) 



(61) 



Summing the first 100 terms in Eqns. (|55)) - (f5"oT) we find from Eq. (|58l) that P = —28A6wo/h 3 so the error in the 
expression for P in Eq. (|5Tj) is less than 5%. The velocity field for £ = h is shown in Fig. [5ja-b). 



The case £/h^>l 



To determine the velocity field in the case f>fcwe again use Eq. (|59|) . Near the centerline parallel to the x-axis, 
i.e. for \y\ < £, we find from Eq. {59) that the wall velocity u(0, 0,z) ~ -AP^f/w 4 so with P = -w^Wq / (4h s ) the 
condition in Eq. (|58"|) is fulfilled. For the transverse velocity u x along the centerline y = we obtain from Eq. ([5"5j) 

■Uz(x,0, z) = Wo/ > -^cos — — . 

n— 1,3,... 



(62) 



This series can be summed 



u x (x,0,z) 



7T 4 h 3 - 6hx 2 + Ax 3 



96 



h 3 



wof- 



(63) 
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The prefactor 7r 4 /96 ~ 1.01, in good agreement with Eq. (|45|) . Similarly, we find from Eq. (|57|) 

u z (x, 0, z) = ^ h ~^ x Wo F (64) 
The prefactor 7r 4 /8 ~ 12.18 in good agreement with Eq. (|46|) . The velocity field for £ 3> /i is shown in Fig. [Sfc-d). 



CONCLUSION 



We have analyzed slow flow in channels with porous walls. A similarity transformation reduces the Navicr-Stokcs 
equations to a set of coupled equations for the velocity potential in two dimensions. We have shown that when the 
Reynolds number and channel aspect ratio is small, an analogy exists between flow in channels with porous walls and 
bending of simply suspended plates under uniform load. If a solution of the problem of a suspended plate has been 
obtained, a problem of viscous motion in a porous tube has thus also been solved. We have applied this result to flow 
in rectangular and triangular channels. Our results provide a general framework for the extension of Berman flow Q 
to three dimensions. 
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